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ON THE ORDER OF MAGNITUDE OF THE COEFFICIENTS IN 
TRIGONOMETRIC INTERPOLATION* 

BY 

DUNHAM JACKSON 

1. Introduction. It is well known that there is a remarkable and far- 
reaching parallelism between the properties of the trigonometric interpolating 
formula determined by the values of a given function at a set of equally spaced 
points, and those of the Fourier's series for the same function, f This paral' 
lelism is closer than the correspondence of sums and integrals in the coef- 
ficients would indicate at first sight. Let /(a;) be the given function, assumed 
to have the period 2ir. Let the Fourier's series, to terms of the pth order, 
be denoted by 

Spix) =^7 + 2^ (on cos nx + 6„ sin nx) . 

^ M=l 

Its coefficients are 

1 r*' 1 f' 

o„ = - I f{x) cosnxdx, bn = - I f{x) sinnxdx. 

'"'Jo IT Jo 

Let the finite trigonometric sum of order p which agrees in value with /(a;) 
at 2p + 1 points equally spaced over the interval from to 2ir be 

(rp{x) = -^ + ^ {Onp cos nx + 6„p sin nx); 
^ 11=1 

in particular, it will be assumed that the points of coincidence are the points 
Xk = 2kir/{ 2p + 1) , k = , 1 , ••• ,2p. The coefficients here are the sums 

2 *E» 2 * 

(1) o„p = zlfixk) cos nxk, b„p = ^ Jlfixk) sin nx*. 

^P -\- i- k-O j:p -\- I J=o 

The resemblance of these to the integrals for a„ and 6„ is apparent; and it is 

* Presented to the Society, December 30, 1919. 

t Cf., e.g., Bocher, Introduction to the theory of Fourier's series, Annals of Mathe- 
matics, vol. 7 (1905-06), pp. 81-152; pp. 132-137; Faber, Ueber sUtige Funktionen, 
Mathematische Annalen, vol. 69 (1910), pp. 372-443; pp. 417-443; D. Jackson, 
On the accuracy of trigonometric interpolation, these Transactions, vol. 14 (1913), pp. 
453-461. 
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322 DUNHAM JACKSON [July 

seen that if n is held fast and p increases without limit, 
lim Onp = On , lim 6„p = 6„ . 

n=ao n=ao 

For any fixed value of p , however, if n is nearly equal to p , the variations of 
cos nx and sin nx are only very imperfectly represented by the values of 
these functions at a set of 2p + 1 points, so that the similarity between Sp 
and (Tp has to be regarded as a consequence of the special properties of the 
trigonometric functions, not merely of the definition of an integral as the 
limit of a sum. 

The present paper is concerned with the order of magnitude of the coef- 
ficients in the interpolating formula, in case the function /(a;) is of limited 
variation, or has a derivative of limited variation. The method of partial 
summation is used to give results analogous to those which Picard* obtains 
for the Fourier's series by the second law of the mean and integration by 
parts. Incidentally, in the case of a derivative of limited variation, a relatively 
simple proof is found for the convergence of the Fourier's series, as well as 
of the interpolating function, to the valuet/(a;) . 

2. Order of magnitude of the coefficients for a function of limited variation. 
Let /(a;) be a function of period 2tr , having limited variation in any finite 
interval, its total variation in the interval from to 27r being V . Then, J 
if On and 6„ are the coefficients in the Fourier's series for/(a;) , as above, 

Kl ^ — , l^nl S — , n. > 0, 

where Ki is independent of n , and may be taken equal to F/tt . This can be 
seen as follows. The function /(a;) — /(O) can be expressed in the form 

/(a;)-/(0) '=fi{x)-f,{x), 

where /i( a;) and/2(a;), the positive and negative variations of f{x) from 
to a:, are non-decreasing functions which vanish for x = 0. Since 

/(2t)=/(0) and fi(2ir)+f2(2ir) = V, 
it follows that 



/i(2x)=/2(2,r) =^. 



The value of o„ is 



a„ = — I /i(a;) cos na;da; — I f2{x) cos nxdx+f(0) I cosnardxj. 

• Traits d'ancdyse, pp. 252-253, 255-256 (second edition). 

t It is immediately apparent from the order of magnitude of the coefficients that the Four- 
ier's series is uniformly convergent, but the proof that the sum of the series is identical with 
/ ( X ) usually requires more attention. 

X Cf. Picard, loc. cit., pp. 252-253. 
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The third integral in the brackets is zero, for n > . By the second law of 
the mean, 



I Ji ( a; ) cos rex ax = ji ( ZTT ) I cos nxdx= -^X = — 



V sin n^ 



where S ^ ^ 27r , and hence 



Similarly, 
and therefore 



I fi(x) cos nxdxl'^ — . 

\Jo \ 2n 





.2w 






./2 ( a; ) cos nx dx 



2/1 ' 



F_ 
nw' 



A corresponding computation for 6„ would seem at first sight to give only 
the inequality 1 6„ | ^ 2F/( nw ) , since the value of /" sin nx dx , extended from 
an unknown lower limit to 27r , might be as great as 2 . For any particular 
value of re, however, the integral defining b„ can equally well be extended 
over an interval of length 27r beginning at the point 7r/(2re) , and so termi- 
nating at a point where cos nx vanishes ; and by expressing f{x) — f{ir/{2n)) 
as the difference of two non-decreasing functions which vanish for a; = 7r/( 2re) , 
it can be shown that 

nw 

An analogous result for the coefficients a„p, hnp is as follows: 
Theorem I. Iff{x) is a function of period 2w which is of limited variation 
in any finite interval, then 

\a„p\ S — , |6„p1 ^— , re >0, 

where Li is independent of n and p . 

If V is the total variation off{x) in the interval from to 2w , the constant Li 
can be taken equal to V . 

It will be convenient to consider the sums 

llfi^k) cos nxk and ^fixk) sin nxk 

of (1) together, as the real part and the coefficient of the pure imaginary part of 
the expression 

2o = Z/(xOe''"'*. 
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Let 

w, = Z e"^** = i. e*- [^■"«^''+»J ( 9 = 0, 1 , . . . , 2p ); 

then 

So =fiXo) ' Wo +f{Zi)(Wi - V)o) +f{Xi){W2 - Wi) 

+ ••• +fiXip) {Wip - W2p_l) 

= [/(a;o) -/(a;i)]wo + [/(ari) -f{xi)]wi 
(2) 

+ ••• + [/(a;2p-i) -/(.T2p)]m!2p_i +/(a;2p) -^25. 

The last term written down vanishes, because* uhp = . Let H he a, number 
such that 

forg = 0, 1, •••,2p - 1. Then 

|So| ^ [\f{xo) -fix^)\ + |/(xi) -f{x,)\ 

+ ••• + |/(a;2p_i) -/(a:2p)i]H, 

or, as the sum in brackets can not exceed V, 

(3) ISol ^ VH. 

It remains to form an estimate of the magnitude of H . 

By applying to Wg the formula for the sum of g + 1 terms of a geometric 
progression, it is seen that 

1 _ g[2(»+l)>nri;(2p+l)] 
(4) Wfl = 1 _ g[2»ir«;Cip+l)] • 

The second term in the numerator is represented by a point on the unit circle 
in the complex plane, so that the absolute value of the numerator can not 
exceed 2. It is evident geometrically (and of course can readily be proved 
analytically) that the absolute value of the denominator is 2 sin [nir/{2p + 1 ) ] . 
Since n/(2p + 1) < 2, and since sin y/y decreases steadily as y increases 
from to 7r/2 , 

nir _ nir sin [mr/{2p + 1)] rnr sin (y/2) 

^^° 2^T1 ~ 2p + r nir/{2p + 1) ^ 2p + l' ir/2 ' 

that is, 

nir 2re 

sm X r—r > 



2p + 1 " 2p + 1 • 



* This is a consequence of well-known properties of the sine and cosine, and, in particular, 
can be inferred immediately from the formula (4) below. 
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It is found consequently that 

1 I < 2p + 1 

and the quantity on the right can be substituted for H in (3). 

Finally, since Onp and fe„p are respectively the real part and the coefficient 
of the pure imaginary part of 22o/(2p + 1) , 

V V 

3. Order of magnitude of the coefficients for a function having a first deriva- 
tive of limited variation. Let it be assumed now that /(a;) has a first deriva- 
tive which is of limited variation in any finite interval, and let V in this case 
stand for the total variation of the derivative in the interval from to 2x. 
By integration by parts,* 

X2lr 1 /»2ir 
f{x)cosnxdx— I f'{x)sinnxdx, 

Xiw 1 /»2ir 

f{x) sin nxdx = - j f (x) cos nxdx, 

the additional terms which would otherwise appear on the right vanishing 
because of the periodicity of the functions involved. Since /'(a;) satisfies 
the hypotheses that were imposed on /(a;) at the beginning of the preceding 
section, 

V 



and 



I /.2T _ V \ C^' 

I f'{x) sin nxdx ^— , I f (x) cos nxdx 



n 



kl^£. |M^£. 



For the case of interpolation, consider once more the sum So , and its repre- 
sentation by the expression (2). It has already been remarked that iHp = 0; 
let the last term in (2), instead of being dropped, be replaced by 

[/(a;2p) -/(^ap+Olwzp, 

with the understanding that Xip+\ = 27r . Then, if the general value of w, 
from (4) is substituted in (2), it appears that 



So = 



L«=0 8=0 J 



1 _ g[2»»'/(?P+<)] 

The first sum on the right reduces to/ (.To) — /(a;2p+.i) , which is zero, because 
• a. Kcard, loc. cit., pp. 255-256. 
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of the periodicity of / ( x ) . In the second sum, let the subscript q be replaced 
by k , and let 

4 = f{xk+i) -f{xk); 

since a;,^! = Xg + 2-kI{2'p + 1 ) , the expression for So reduces to 

.[2n.ri/(2p+I)] 2p 



(5) 



2o 



1 



H)l Z^"* 



.7. 

[2»«i/c2p+i)] Z^"* e 
e *=0 



This formula corresponds to the result of the integration by parts in the case 
of the Fourier's coefficients. 

In the factor before the sign of summation on the right, the absolute value 
of the numerator is 1 . The denominator is the one that was considered in the 
preceding section; its absolute value is 2 sin [n7r/(2p + 1)], and is greater 
than 4n/(2p + 1) • Consequently 



|So| ^ 
By the mean value theorem, 



Erf, 



k e" 



2-n- 



dk = (xk+i - Xk)f'{^k) = 2jr_n/'(^it) ' 3:k < ^k < Xk+1. 

Let a function (^ ( x ) be defined so as to be equal to /' ( f t ) f or a; = % , Jfc = , 
1, ■■• , 2p, equal to/'(^o) for x = 2w , linear between any two successive 
points 2kTr/{2p + 1), and periodic with the period 2ir. Then <t>{x) is of 
limited variation in any finite interval, its graph being made up of a finite 
number of linear segments, and so satisfies the conditions imposed on/(a;) 
at the beginning of the preceding section. Its total variation from x = 
to X = 27r is 

lf(^l)-/'(?0)| + |/'(^2)-f($l)| + --- 

+ If (?2p) -/'(f2p-l)| + If (^0 + 27r) -f (f2p)|, 

and this quantity, while not necessarily equal to F , is at any rate not greater 
than V. Hence the reasoning of the earlier section is applicable, and it 
appears that 



Z<A(^Oe"''' 



^ V 



2p + l 

2n ' 



From this can be deduced successively the inequalities 



Hdi 



k e" 



a. 



n ' 



So ^ 



{2p + l)Vw 



4/1^ 



Ftt 

K\ =2^2- 
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The conclusion may be stated as follows :* 

Theorem II. If f{x) is a function of period 2ir which has a first derivative 
f {x) of limited variation in any finite interval, then 

I Unp 1=2' I "np 1=2' II' ^ yJ , 

where L2 is independent of n and p . 

If V is the total variation off {x) in the interval from Q to2ir , the constant Li 
can he taken equal to Vt/2 . 

4. Proof of convergence. If /(.t) satisfies the hypotheses of the preceding 
section, the inequalities found there for the coefficients a„ , bn make it clear 
that the Fourier's series for /(a;) is uniformly convergent. It is not so im- 
mediately obvious from first principles that the sum of the series must be/(a;) . 
On the other hand, since a„p and bnp depend on p as well as on re, it is not 
evident, even as a result of Theorem II, that the interpolating function cTp{x) 
approaches a definite limit as p becomes infinite; but it is readily seen that 
if (Tpix) does converge uniformly to a limiting function g{x), then g{x) 
must be identical with /(a;) . More explicitly, ii <Tp{x) converges to the limit 
g{x) when p becomes infinite through all integral values, it will converge to 
the same limit if p is restricted to the values (3'" — 1 )/2 , so that 2p + 1 = 3™ , 
m = 1,2,3, • • • . But in the latter case, if x is any particular value of the 
form 2A:7r/(3'*) , where 11 is any positive integer and k is any integer whatever, 
(Tpix) will be equal to f{x) for m = fi and for all succeeding values of m, 
because of the properties oi ffp{x) as an interpolating function, and therefore 
the value of g{x) at this point must be the same as that ot f{x) . That is, 
g{x) is equal to f{x) at all points x = 2A;ir/(3''); and as both /(.r) and 
g{x) are continuous, and these points are everywhere dense, g{x) must be 
equal tof{x) identically. It will be shown that <Tp{x) converges uniformly 
to the same limit as Sp{x) , and then it will follow, by what has just been 
said, that both approximating functions converge! to the value /(a;) . 

Let 

g{x) = Umspix); 

p=00 

.also, let (Tnpix) be the sum of the first terms of ap{x), through those of 

* In connection with the value of L, , cf. § 5 below. 

t The truth of this conclusion is well known, for the interpolating function as well as for 
the Fourier's series, with more general hypotheses as to the function fix); cf . Faber, loc. cit., 
and also de la Vallee Poussin, Sur la convergence des formules d'interpolation enire ordonnees 
equidislanles, Bulletins de I'Academie royale de Belgique, Classe 
des Sciences, 1908, pp. 319^03. The purpose of the present section is to answer 
the natural question what connection there is between the order of magnitude of the coef- 
ficients and the convergence of the interpolating function, and incidentally to call attention 
to the simplicity of the resulting convergence proof. 
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order n , 



o'np (a;) = -^ + 23 (ffltp cos kx + bkp sin kx) . 



The convergence of Sp{x), as already noted, follows from the inequalities 
derived in the preceding section for a„ and fe„ . In consequence of this con- 
vergence, and of the inequalities for o„p and Kp, there will exist for every 
positive € a subscript n such that 

\gix) - s„ix)\ <g, \<Tp{x) -cr„p{x)\ <g, 

both inequalities holding for all values of x , and the latter one for all values 
of p > n . If this value of n is held fast, and p is allowed to increase without 
limit, akp and bkp will approach Uk and bk respectively, for all values of k ^ n, 
(fupix) will uniformly approach 5„ (a;) , and, when p is suflSciently large, 

\Sn{x) - <T„p(x)\ < 2 

for all values of x . By combination of the last three inequalities, it is seen that 

\gix) - (Tp{x)\ < e 

as soon as p is sufficiently large, and this completes the proof of convergence. 
5. Alternative method for a function having a derivative of limited variation. 
If the convergence of the Fourier's series to the value /(.r) , in the case that 
f{x) has a first derivative of limited variation, be assumed as known at the 
start, inequalities for the order of magnitude of the coefficients in the inter- 
polating formula can be obtained without the use of partial summation. 
When/(.rt) in (1) is replaced by its expression in terms of the Fourier's series, 
it is found that* 

fflnp = a» + 2^ ( C'\(.2p+l)+n + aA(2p+l)-»i ) , 

inp — K + Zl (^X(2p+1)-H» — bx(2p+l)-n) • 
A=l 

On the other hand, by the inequalities at the beginmng of § 3, 

I I <:^ \h \ ^ — 

where K2 is a constant. Furthermore, since n ^ p, 

X(22J + 1) +n > (2X + l)re, X(22) + 1) - n > (2X - l)n, 
* Cf., e.g., de la Vallee Poussin, loc. cit., p. 383. 
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1 a^(2p+i)+n I = ( 2X + 1 )2 n^ ' 1 "^(2^+1)+" I = {2\ + \Yn" 

1 aM2p+l)-„ I = ( 2X _ 1 )2 ;i2 ' 1 ''A(2p+i)-n 1 = ( 2X - 1 )2 n^ ' 



so 



that 



ttn 



: „2 \^i +32 + 52+ •••^ - ,^2. |0np|=^2. 



L'2 being independent of n . 
Since* 

1+1+1+... =1! 
■1 + 32 + 52 + 3 » 

and since, by the preceding section, K2 may be taken equal to V/ir , the value 
F7r/4 is found for L'2 , which is an improvement over the value of L2 given in 
the statement of Theorem II. 

li f{x) possesses a derivative of order higher than the first having limited 
variation, the first derivative is certainly of limited variation, and it follows 
from § 4, even if it is not assumed otherwise, that the Fourier's series con- 
verges to the value f{x). The method just outlined can be used, with the 
stronger inequalities for a„ and 6„ which are obtained by integration by parts 
in this case, to yield corresponding inequalities for Onp and 6„p . 

6. The method of partial summation for a function having a second deriva- 
tive of limited variation. Although the existence of higher derivatives can 
be taken into account in the manner indicated at the end of the preceding 
section, it is not without interest to observe that the method of §§ 2, 3 can be 
made to apply to this case also. The extension is not entirely trivial, but it 
will perhaps be sufficient to carry it one stage further. Let f" {x) , then, be 
everywhere defined, and of limited variation in any finite interval. It is to 
be shown that 

I ^np 1 = -^ , 1 Onp 1 = "^ , 

where is is independent of n and p. The more exact determination of the 
constant is of secondary interest. 



Let it be understood now that 



2kT 

Xk = 



2p + l' 

for all integral values of k , whether contained in the interval from to 2p + 1 
or not; let 

i'k =fixk+i) -f{xk) = dk, 



* Cf ., e.g., Bromwich, Theory of infinile series, p. 187. 
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an accent being added to the symbol dk previously used, to indicate that these 
are first differences, as distinguished from the second differences now to be 
introduced; and let 

du = dk+i — du . 

The coefficients a„p and 6„p , apart from the factor 2/( 2^ + 1 ) , are the com- 
ponents of the complex quantity So already considered. If 

k=0 

it follows from (5) that 

g(2»ir«/(2p+l)] 
2o = J _ gr2r.,.7(2p+l)] 2l . 

The sum Si can be rearranged by partial summation, after the same manner 
as So . It is found that 

Si = d'oWo + d'l (wi — Wo) + • • • + d2p{w2p — W2p-i) 

= {d'o — d'i)wQ + {d\ — d2)wi + ■• • + ( d'lp-i — d'lp ) w^p-i + d'lp Wip-, 

since wip = , the last term may be replaced by ( d'^p — d'ip+i ) W2p . In 
consequence of (4), furthermore, 



J r 2p 2p -| 

w«2p?T7] T,id',- d;+i ) -Hid',- d',+i ) e"^^'* . 

L?=0 9-0 J 



'"^ ~ I _ g[2niri/(2p+l)] 

The first sum on the right reduces to d'o — rfjp+i, and this is equal to zero, 
because rf2p+i is the same as d'o . Hence 



t2nff«/(2p+I)] 2p 



Si — , ^[2»ir(/(2p+nl 2^ "'k ' 



,"»•** 



1 ^[2»ir(A2p+l)l Z^ "'k B 

and 

/ g[2n,r</(2p+l)] \ 2 2p 
2^0 = \\ — g[2^>'</(2P+l)] ) ^ * ^"^* ■ 

It is clear that a corresponding relation can be obtained for differences of any 
order. By reason of the inequality previously used for the factor in paren- 
theses, 



<,ni« 



It remains to discuss the order of magnitude of the sum on the right. 

By one more partial summation, which it is not necessary to give in detail, 
it is seen that 



2P 

E d'k e" 

t=0 



^(|rfr-4'l + \d!^ -d'^\ + ••• + \d'ip - d'ip.,\) H , 
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where H , as in § 2, is an upper bound for | w, | , and may be taken equal to 
(2p + l)/(2n) . So, if X stands for the quantity in parentheses, 

(6, 1^.1^^^!^'^. 

and the problem is further reduced to that of treating the expression X . 

If a function /(x) , possessing a second derivative throughout an interval 
a '^ X ^ c, vanishes for a; = a and for x = c, and also at an intermediate 
point X = b, the first derivative /'(a:) must vanish at some point between 
a and b, and again between b and c, and so f" (x) must vanish somewhere 
between a and c . If / ( a; ) is not required to vanish at the points a,b , and c , 
but A -{■ Bx + Cs? is the polynomial of the second degree agreeing in value 
with /(a;) at these three points, the second derivative of 

fix) -{A + Bx + C3?) 

must vanish somewhere between a and c , or, in other words, there is a point 
at which /" ( a; ) = 2C . Now let it be assumed further that b is midway between 
a and c, so that b = a + h, c = a -\- 2H, say. Then the second difference 
of /(a;) for the three points, 

j{x + 2h) -2fix + h)+fix), 

is the same as the corresponding second difference formed for the function 
A + Bx + Cx? , and the latter second difference is equal to 2Ch^ . That is, 
since 2C =/"({) , where | has some value between o and c , the second differ- 
ence of /(a;) is equal to* h^f" {^) . 

In the problem under discussion, where h — 2x7 (2p + 1) , 

where fit is some point between x* and Xk+i . Consequently 
^ = (-2~^2[i/"(?i) -/"(«o)| + |/"(f2) -/"($i)| 

+ --- + i/"(«2p)-r(i2p-i)|]. 

if it could be said that ?o < f i < ^2 < • • • , the quantity in brackets could 
be brought into immediate connection with the total variation of /"(a-); 
but it is not clear that the |'s do occur exactly in the order of their subscripts. 
However, since .To < ?o < a;2 , Xi < |i < xs , etc., there is no doubt that the 
pair of points ( fo , ?i ) precedes the pair ( {3 , ^4) , etc., and as the order of the 
terms in any one difference \f" ( ^k+i ) — /" (I* ) | is immaterial, it can be 
* Cf. Markoff, DifferenzeriTechnung (Leipzig, Teubner, 1896), p. 21. 
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asserted that for present purposes the points ^o, ^i, f 3 , ?4 , ^e , ^7 , • ■ • , are 
essentially in order of algebraic magnitude. Moreover, while there may be a 
question whether the points ^0 > • • • . ^2p are all included in an interval of 
length 2ir , they are at any rate all included in an interval of length iir . So, 
if V is the total variation of/" (a;) in an interval of length 2-w , it is certain that 
If'Ui) -/"(?o)| + |/"(?4) -r'{b)\ + If'C^r) -/"(^6)| + ••• ^2F. 
Similarly, 

|/"(?2)-/"(^l)| + l/"($5)-/"(|4)| + ---S2F. 

i/"(^3) -r(^2)i + ircie) -/"(?6)i + ••• ^2F. 

It would be possible to do something toward replacing these by a stricter set 
of inequalities, but that point will not be pressed. By combination of the 
three as they stand, 

l/"(^i) -/"(?o)l + l/"(^2) -/"(^i)| + ••• ^6F, 

247r'r 
^ = (2p + l)2- 
With (6), this gives 

Btt^F 2p + l 
|2o| =-4 ^^. 

Finally, because of the factor 2/(2p + 1) in Unp and 6„p, 



dnii = 



2^3 , li/npl = 2^3 



the factor Stt^ F/2 being independent of n and j> . 

There is no further essential difficulty in carrying the discussion on to 
derivatives of higher order. 

The University op Minnesota, 
Minneapolis, Minn. 



